Abstract 



A search has been made for the decays B + — » 4>ir + and B — * 4>ir 1 in a 
data sample of approximately 232 xlO 6 BB pairs recorded at the T(45) 
resonance with the BaBar detector at the PEP-II B-meson Factory at 
SLAC. No significant signals have been observed, and therefore upper 
limits have been set on the branching fractions: B(B + — > 4>tt + ) < 2.4 x 
1(T 7 and B(B° -» (pn ) < 2.8 x 10~ 7 at 90% probability. 
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A search has been made for the decays B + — > (f>ir + and B — > 6n° in a data sample of approxi- 
mately 232 xlO 6 BB pairs recorded at the T(4S) resonance with the BABAR, detector at the PEP-II 
B-meson Factory at SLAC. No significant signals have been observed, and therefore upper limits 
have been set on the branching fractions: B(B+ -> 0tt+) < 2.4 x 10~ r and B(B° -> fin ) < 2.8 x 10~ r 
at 90% probability. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

The measurements of B — > <j)K and B — > 07r decay 
rates are important because they are sensitive to con- 
tributions beyond the Standard Model (SM). In partic- 
ular, the latter is strongly suppressed in the SM, and a 
measurement of B(B — > </>7r) > 10~ 7 would be evidence 
for new physics, for example supersymmetric contribu- 
tions The study of the processes B + — > <f>iT + [2( and 
B° — > <fiir° is also important to understand the theoret- 
ical uncertainties associated with measurements of CP 
asymmetries in B° — > <fiK° decays. The B — > <f>ir decay 
amplitudes are related to the sub-leading terms of the 
B° — > (f>K° decay amplitude Q and can therefore pro- 
vide stringent bounds on possible contributions to the 
time-dependent CP asymmetry in B° — ► <pK° p|, an- 
other probe of new physics effects in B decays. 

In Fig. ^ w e show the leading order Feynman diagram 
for the B — > </>7r decay and a sub-leading diagram for 
B — > (f>K decay. 




4> 
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FIG. 1: Feynman diagram for B — > <ji>7r and B — > <j>K. 

Previous searches for these decay modes have been re- 
ported by BABAR and CLEO H H The results pre- 
sented here are based on data collected with the BABAR 
detector || at the PEP-II asymmetric-energy e + e _ col- 
lider located at the Stanford Linear Accelerator Cen- 
ter. An integrated luminosity of 211 fb _1 , corresponding 
to (231.8±2.6) x 10 6 BB pairs, was recorded at the T(AS) 
resonance (center-of-mass energy ^fs = 10.58 GeV). 

Charged particles from the e + e~ interactions are de- 
tected, and their momenta measured, by a combination 
of five layers of double-sided silicon microstrip detectors 



(SVT) and a 40-layer drift chamber (DCH), both op- 
crating in the 1.5-T magnetic field of a superconduct- 
ing solenoid. Photons and electrons are identified with 
a CsI(Tl) electromagnetic calorimeter (EMC). Further 
charged particle identification (PID) is provided by the 
average energy loss (dE/dx) in the tracking devices and 
by an internally reflecting ring imaging Cherenkov detec- 
tor (DIRC) covering the central region. 

Monte Carlo simulation is used to evaluate background 
contamination and selection efficiency. Signal and back- 
ground Monte Carlo samples are generated with EvtGen 
llfj. The detector response is simulated with GEANT4 
[ll| and all simulated events are reconstructed in the 
same manner as data. 

We reconstruct B meson candidates through the de- 
cays 4>tt + or (f>ir°, with <f> — > K + K~ and 7r° — > 77. All 
kaon candidate tracks in the reconstructed decay chains 
must satisfy a set of loose kaon identification criteria 
based on the response of the DIRC and the dE/dx mea- 
surements in the DCH and SVT. In both decay modes, 
all the tracks coming from the fully reconstructed B are 
required to originate from the interaction point. A pair 
of oppositely-charged kaon candidates is considered as a 
4> candidate if its invariant mass is within 15MeV/c 2 of 
the nominal <f> mass value (1019.5 MeV/c 2 [l2()- This is 
about three times the observed width in the K + K~ in- 
variant mass spectrum. A pair of energy deposits in the 
EMC, each of which is isolated from any charged track 
and has the lateral shower shape expected for photons, 
is considered as a tt° candidate if both the deposits ex- 
ceed 40 MeV in the laboratory frame and the associated 
invariant mass of the pair is between 110 MeV/ c 2 and 
160 MeV/ c 2 (about three times the observed width in the 
77 invariant mass spectrum). B meson candidates are 
made by combining <fi candidates with a charged track 
or a 7T° candidate. We do not apply any particle identi- 
fication criteria on the track which comes directly from 
B meson decay (primary track) at this stage, so for the 
charged mode we reconstruct B + — > <f)h + (h + = ir, K) 
events. This allows us to study the B + —> (j)K + signal, 
which is the largest background coming from B decays. 

Two kinematic variables are used to discriminate be- 
tween signal B decays and combinatorial background: 
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the invariant mass of the reconstructed B meson can- 
didate, m B and m m i ss = \J (q e + e - ~ 9b) 2 , where q e + e - 
is the four momentum of the initial e + e~ system and 
qb is the mass-constrained four momentum of the re- 
constructed B meson candidate. By construction, the 
linear correlation between m m i SS and tub vanishes. Com- 
pared to the kinematic variables AE = E B — \\J~s and 

mES = \J \ s ~ P*b ( wn ere s = q\+ e - and the asterisk 
denotes the e + e _ rest frame), which were used in the 
previous BABAR analysis of these modes 0, the present 
combination of variables has less correlation and bet- 
ter background suppression. The distribution of m_B 
peaks at the nominal B mass value |l2j . with a width 
of about 20MeV/c 2 for 0tt+, and about 40MeV/c 2 for 
(j>Tr° , with a low-side tail due to energy leakage from 
the EMC. The resolution on m miss is about 5MeV/c 2 , 
dominated by the beam-energy spread. We require tub 
to be within 150 MeV/c 2 of the nominal B mass and 
5.11 GeV/c 2 < m miss < 5.31 GeV/c 2 . The region 
Wmiss < 5.2 GeV/c 2 is used for background character- 
ization. 

The dominant background comes from combinatorial 
e + e~ — > qq (q = u, d, s, c) continuum events. They tend 
to be jet-like in the center-of-mass (CM) frame, while B 
decays tend to be spherical. To exploit this characteristic 
for discriminating against continuum background, we use 
the ratio L2/L0, where Li is defined as 



Li 



Eh 



(i) 



where pk is the momentum of particle k, and 9k is the an- 
gle between pk and the thrust axis of the reconstructed B 
meson evaluated in the CM frame. The sum runs over the 
charged and neutral particles of the event not assigned 
to the B meson. We require L2/L0 < 0.55, which sup- 
presses the continuum background by more than a factor 
of 3, while retaining about 90% of the signal. We require 
I cos 9* B I < 0.9, where 9* B is the angle between the B 
candidate momentum and the e + momentum in the CM 
frame. For B candidates the probability density function 
of 9* B is proportional to sin 2 9* B , whereas for continuum 
events it is nearly uniform after acceptance. We select 
events for which one B is reconstructed as B + — > <fih + 
or B° — > 4>tt° and the other B is only partially recon- 
structed [13. We define At to be the difference between 
the proper decay times of the B mesons and a At the un- 
certainty associated with it. We require, in the case of 
B° -> 07T° only, |A*| < 20 ps and a A t < 2.5 ps. These 
requirements on At and aAt retain about 92% of the sig- 
nal, while removing about 15% of the continuum events. 
The r.m.s. At resolution is 1.1 ps for the events that 
satisfy these requirements. After the application of these 
selection criteria on Monte Carlo simulated events, the 
efficiencies for 07r + and </>7r° signal are (37.1 ±0.1)% and 
(29.5 ± 0.8)% respectively. The average candidate mul- 



tiplicity in events with at least one candidate is ~ 1.005 
for both decay modes. If more than one B candidate 
is reconstructed in an event, we choose the one with the 
<t> — > K + K~ invariant mass closest to the nominal tf> mass 
value ,12], for B + — > (f>ir + decays. For B° — > (ftir decays, 
we choose the candidate with the 7r° — ► 77 invariant mass 



t_L7 

closest to the nominal 7r u mass value |l2|. These criteria 
produce no bias in the shape of the other event variables 
used in the maximum likelihood fit described below. We 
select 10990 and 2732 events in the <fih + and ipir analyses 
respectively. 

A possible background to the <f> — > K + K~ decays 
comes from the S-wave production of the K + K~ system 
(B — > (K + if~)s-wave7'" decays) with contributions com- 
ing predominantly from resonances such as / (980) and 
a (980). Using samples of simulated decays of B mesons 
equivalent to nearly five times the size of the data sample, 
we found that all the other B decay modes give negligible 
sources of background. To discriminate against S-wave 
background in the maximum likelihood fit, we use the 
helicity of the K + K~ system, in terms of the cosine of 
the angle 9h between the K + candidate and the par- 
ent B meson flight direction in the K + K~ rest frame. 
The helicity probability density function is proportional 
to cos 2 9h for the signal, and is uniformly distributed for 
the S-wave background. Further discrimination is pro- 
vided by the K + K~ invariant mass distribution, uikk, 
which peaks at the <f> mass for the signal, while it peaks 
at lower values for the S-wave background. 

In the case of charged B decays, we exploit the 
Cherenkov angle 9 C measured in the DIRC for the pri- 
mary track, in order to determine simultaneously the 
yields of B + — ► cf>ir + and B + — > c/)K + decays and the 
yields of the two corresponding B + — > (K + K ~)s~wnveh + 
(h = 7r, K) background components. 

Signal and background yields Ni, where i denotes sig- 
nal, continuum, and S-wave background, are extracted 
using an extended maximum likelihood fit with the like- 
lihood function: 



C 



exp 



n 



^NiPiixjiSi) 



(2) 



where N is the total number of events entering the fit. 
The probabilities Vi are products of Probability Den- 
sity Functions (PDF) for each of the independent vari- 
ables x = {m m i SS , tub, L2/L0, rriKK, cos 9h}- In the case 
of B + — * cj)h + the variable 9 C is also used in the fit. 
The di are the parameters of the PDFs for x. The con- 
tinuum parameters are allowed to vary, except for the 
m m i ss end-point. All other parameters cti are fixed to 
their values derived from data control samples. These 
are varied within their uncertainties to evaluate the sys- 
tematic error. By minimizing the quantity — ln£ in two 
separate fits, we determine the yields for <j)ir + and tfiir . 



6 



There are three B backgrounds to B + — > 4>ir + decay 
(B+ -> {K+K-) S - wavc h+ and B+ -> <j>K+), while only 
B° -> (i^+i^-) s _ wavc 7r° contributes to the B° -> </>7r° 
mode. All the yields in Eq. [3 are allowed to fluctuate to 
negative values in the fits. 

The distributions of L? I Ln and cos Or are described by 
a parametric step function [14j and a second-order poly- 
nomial respectively. We use a Gaussian for the m miss 
distribution for cj)ir + signal and S-wave components, a 
Gaussian with exponential tails for the tub distribution 
for 07r + signal and S-wave components, and for both 



signal and S-wave components. 
'hubs distribution we use the function 



m m i ss and tub for <pn" 
For the continuum m n 
x\J 1 — x 2 exp [— £(1 — a; 2 )] , with x = 2m m i SS /v / i and £ a 
floating parameter. The ttikk invariant mass distribu- 
tion is described by a relativistic Breit-Wigner function 
for signal, a relativistic Breit-Wigner plus ex pon ential for 
the continuum background and a Flatte jlfii function 
for the S-wave background. The Flatte function takes 
into account the coupling of the scalar resonances to the 
7r + 7r~ and K + K~ channels 17]. 

The Cherenkov angle C PDFs are obtained from a 
large data sample of D*+ -> D n + (D° -> A'-tt+) 
decays where A^/tt* tracks are identified through the 
charge correlation with the 7r ± from the D*^ decay. The 
PDFs are constructed separately for A + , K~, ir + and 
7r~ tracks as a function of momentum and polar angle 
using the expected values of 9 C , and its uncertainty. 

Using a large number of simulated experiments, we 
find that the usual maximum likelihood fitting technique 
does not provide an unbiased estimate of the true values 
of signal and S-wave yields (As and As- waV c) because of 
the non-Gaussian shape of the likelihood function when 
the yield is very small. Therefore we use a Bayesian sta- 
tistical approach to obtain a modified likelihood function 
L(N S ): 

L(N S ) = N dA s _ wavc £(A s ,A s _ wavo ), (3) 
Jo 

where the normalization Ao is such that 
J dNsL(Ns) = 1. The two dimensional likeli- 
hood C(Ng, As-wavc) is given at each point on the 
As- As -wave plane by the function defined in Eq. 
maximized with respect to all of the other fit variables. 
When seeking the central value for the branching fraction 
we take the median of L, with the lower limit replaced 
by — oo and As unrestricted. This is because we find 
from simulations that in the case of very low yields, the 
median provides a less biased estimator of the true value 
of As than the maximum of L. We correct the central 
value of the branching fractions for the residual biases. 
When calculating upper limits, we impose the a priori 
constraints As > and As- W avc > 0. 

Fig. |21 shows the m m i ss and vtlb distributions for data, 
with the PDF corresponding to the maximum likelihood 



fit overlaid. We do not observe evidence for either B + 
(f>ir + or B° — > <fiir° decays. 



TABLE I: Signal yield (evaluated as the median of the like- 
lihood), detection efficiency e (the uncertainty includes both 
statistical and systematic effects), measured branching frac- 
tion B with statistical error, after the correction for the fit 
bias has been applied, for the two decay modes considered 
and upper limit at 90% probability. 





B + -> <^7T + 


do j. o 
a — > <pn 


Yield 


-1.5 ±5.9 


4.0 ±3.5 


e(%) 


37.1 ±0.1 


29.5 ±0.8 


B(w~ 6 ) 


-0.04 ±0.17 


0.12 ±0.13 


UL(23)(l(r 7 ) 


2.4 


2.8 



The signal yields, extracted from the median of the 
likelihood L(N S ) (Eq. 0), are reported in Table [Q In 
the case of B + — > <fih + , we also measure the N,p K + yield, 
which is found to be compatible with the expectation 
from published branching fractions 

The branching fraction B is calculated from the ob- 
served number of signal events as 



A s 



e-N BB 



B(<t>^K+K~) 



(4) 



where N B § is the number of BB pairs produced and 
e is the reconstruction efficiency for the B candidates. 
In Eq. 0| we assume equal branching fractions for Y(4S) 
decays to charged and neutral B-mcson pairs 18]. 

The systematic uncertainties are summarized in Ta- 
ble[n] The uncertainty arising from the lack of knowledge 
of continuum background PDFs is part of the statistical 
error since the background parameters are free to vary 
in the fit. The uncertainty on the signal PDFs repre- 
sents the dominant error. We estimate it by using sim- 
ulated and high-statistics data control samples of B + — > 

) 



r + D o ( D o ^ K+TT-) and B° -> tt+ D~ (D~ 



events. In order to estimate the systematic uncertainty 
on tub for B° — * (fm we use a data control sample of 



TABLE II: Summary of systematic uncertainties contributing 
to the total error for the upper limit on the branching fraction. 
They are given in units of 10 -8 . 



B° 



PDF Uncertainty 


+ 1.9 
-2.8 


+3.6 
-4.2 


PID Efficiency 


0.1 


0.1 


Tracking Efficiency 


0.1 


0.2 


7T° Efficiency 




0.1 


L 2 /L Cut 


0.1 


0.3 


BB Pair Counting 


0.1 


0.2 


Interference Effects 


0.3 


0.6 


B(cj>^K + K-), B(ty° ^77) 


0.1 


0.1 


Total 


+2.8 
-3.6 


+3.7 
-4.3 




5.35 5.4 , 

m B (GeV/c) 



5.35 5.4 , 

m B (GeV/c ) 



5.35 5.4 , 

m B (GeV/c) 



FIG. 2: Distribution of m m i ss (top) and mj (bottom) for reconstructed B° — » 4>ir° (left), B + — > 4>tt + (middle) and B + — > 4>h + 
(h = {7r,if}) (right), after applying a requirement on the ratio of signal likelihood to signal- plus-background likelihood to 
enhance the signal. The curves are projections from the likelihood fit for total yield (solid line) , for the continuum background 
(fine dashed line) and for continuum plus S-wave component (dashed line). For B + — > <f>h + decay we do not apply any particle 
identification criteria and assign the pion mass to the primary track. For this reason, the tub distribution for B + — > (j>K + 
events is shifted with respect to the nominal B + mass (positive bump in bottom middle and bottom right plot), while it peaks 
at the nominal value for B + — > 4>ir + events (negative bump in bottom middle plot). 



B + — > h + ir° events. The control channels have event 
topologies similar to those of B + — > <ph + and B° — ► 07r°. 
We use them to determine the signal PDF parameters 
and take the difference in yields found by varying these 
parameters within one standard deviation as the system- 
atic error. The second most important error comes from 
the uncertainty on the efficiency s. The track detection 
efficiency uncertainty is estimated to be 0.8% per track 
from a study of a variety of control samples, such as 
t — * 3-track decays. We assign 0.5% uncertainty on the 
kaon identification efficiency. The uncertainty on the re- 
construction efficiency for the ir° is 3%, as measured in a 



large sample of r 



7r 7T decays coming 



P »t, P 

from e + e~ — > t + t~ . We assign a 1.8% uncertainty on the 
L2 1 Lq cut efficiency, estimated by the difference between 
Monte Carlo and data control samples, 1.1% on the total 
number of T(4S) — > BB decays in the sample and 1.2% 
on the knowledge of B(n° — ► 77) and B(<j> — > K + K~). 
We estimate the systematic error introduced by the ap- 
proximation of ignoring interference effects between the 
<f> and the K + K~ S-wave components by varying the rel- 
ative strong phases and taking the largest observed vari- 
ation as the error. In this study we include the /o(980) 
resonance and a non-resonant component, whose contri- 
bution is taken from a B + — > K + K~K + Dalitz plot 
measurement by the Belle Collaboration^^. The result- 
ing uncertainty is 4.4% for both modes. 



Under the assumption that N B g and e are distributed 
as Gaussians, we obtain a likelihood function, Lg, for 
the branching fraction, B, based on by convolving 

the likelihood (L in Eq.|3) with the distributions of N B § 
and e. We also include the additional uncertainty coming 
from the systematic error on the signal yield. The result- 
ing likelihood is shown in Fig.Olfor each of the two decay 
modes. In the plots, the upper boundary of the dark re- 
gion represents the 90% probability Bayesian upper limit 
Z?tjl, defined as: 



L B {B)dB 



9 

10 



L B (B)dB 



(5) 



We determine B{B+ -> 0tt+) < 2.4 x 10~ 7 and B(B° -> 
cfm a ) < 2.8 x 10~ 7 . 

We compute the central values for the branching frac- 
tions by correcting the fitted signal yields for the fit 
bias, estimated using a large number of simulated ex- 
periments, and including a systematic uncertainty equiv- 
alent to half the fit bias. This error corresponds to a 
shift of —0.8 and —0.4 events in the signal yield, and 
— 1.9 x 10~ 8 and —1.2 x 10~ 8 in the branching fraction 
for B + — > 07r + and B Q — > </>7r°, respectively. Without 
taking into account the a priori knowledge of Ns > 
and -ZVs-wave > 0, and integrating the likelihood in 
Eq. [3] around the median, we obtain as 68%-probability 
regions B{B+ -> cjm + ) = (-0.04 ±0.17) x 10~ 6 and 
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BF ( <twc + )[li BF(<|m )[10" 6 ] 



FIG. 3: Likelihood distribution, L B (B), for B(B + -> 0tt+) 
(left) and B(B° — > 4>tt°) (right) in arbitrary units. The upper 
boundary of the dark region represents the 90% probability 
upper limit. 



B(B° -> </>7T°) = (0.12 ±0.13) x 10~ 6 . The results are 
summarized in Table [i] 

In summary we have searched for B + — > </>7r + and 
B° — > 07T° decays in a sample of 232 million BB me- 
son pairs. We find no evidence of signal and therefore 
we place upper limits B{B + — ► cf>ir + ) < 2.4 x 10~ 7 and 
B(B° -> </>7T°) < 2.8 x 10~ 7 at 90% probability These 
limits are more stringent than earlier results^, 0, an d 
they supersede our previous publications 13, U They are 
consistent with existing SM predictions |lj. 
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